How colloidal suspensions that age are rejuvenated by strain application? 
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We present here a microscopic study of the effect of shear on a dense purely repulsive colloidal 
suspension. We use Multispeckle Diffusing Wave Spectroscopy to monitor the transient motions of 
colloidal particles after being submitted to an oscillatory strain. This technique proves efficient to 
record the time evolution of the relaxation times distribution. After a high oscillatory shear, we show 
that this distribution displays a full aging behavior. Oppositely, when a moderate shear is applied 
the distribution is modified in a non trivial way. Whereas high shear is able to erase all the sample 
history and rejuvenate it, a moderate shear helps it to age. We call this phenomena overaging. 
We demonstrate that overaging can be understood if the complete shape of the relaxation time 
distribution is taken into account. We finally report how the Soft Glassy Rheology model accounts 
for this effect. 

PACS numbers: 07.60.-j, 78.35.+C, 81.40.Cd 



I. INTRODUCTION 



The study of dense colloidal suspensions still raises 
many open questions. One interesting feature is that 
some of theses systems remain out of equilibrium for any 
accessible experimental time. For instance, if the pair po- 
tential of interaction has an attractive part, the particles 
can aggregate in an always evolving physical network [Q . 
If the potential is purely repulsive, an increase of concen- 
tration leads to a jamming transition and the suspension 
can reach an amorphous state or glassy phase More 
precisely, in this last case the characteristic relaxation 
time of the system increases dramatically around a crit- 
ical value of the volume fraction ^c- For large enough 
concentrations, the relaxation processes rapidly become 
longer than any experimental time. The system is then 
macroscopically pasty and can display a non-stationary 
behavior. A direct consequence of the extremely long 
relaxation times is that all the past history of the sam- 
ple has to be taken into account. A key trick to obtain 
reproducible results is to perform experiments on a sam- 
ple with the exact same history. As far as rheological 
measurements are concerned, it is known for long that a 
high preshear can erase all past memories in many sys- 
tems. It is thus often used as a trick to reset the sample 
history provided that it does not damage it irreversibly. 
This procedure is reminiscent of the process of thermal 
quench used to erase the history for other structural or 
spin glasses [|[ ^. This similarity let envision a com- 
parable role played by shear and temperature on a mi- 
croscopic level. It has recently been shown on different 
'pasty' colloidal suspensions that their macroscopic be- 
havior after a rheological quench have the same qualita- 
tive feature than other glasses after a temperature quench 
[H, Ij. However, measurements of the effect of shear on 
the dynamic at a microscopic level are still lacking. In 
this paper we study a colloidal glass with Multispeckle 
Diffusing Wave Spectroscopy (MSDWS), monitoring the 
motions of the colloidal particles after that the sample 
underwent various strain histories. We will show that 



the mechanical perturbation acts in a dual fashion on 
the microscopic dynamic. All the observed effects will be 
discussed within the Soft Glassy Rheology (SGR) model 



II. BACKGROUND 

In this section we will recall some general results on 
glasses. We will mainly focus on results obtained on 
spin and structural glasses upon a change of tempera- 
ture. However, most of theses results are similar to that 
observed in colloidal glasses if a temperature decrease 
is substituted by an increase of concentration ||, 
Spin and molecular glasses display a qualitative change in 
their microscopic dynamic when the temperature is low- 
ered around the glass transition temperature. Actually, 
the distribution of relaxation times splits into two dis- 
tinct families when the temperature is decreased (see eg 
on the one hand, short distance motions and vi- 
brations of the particles can be describe by fast individual 
modes called P modes. Their sensitivity to temperature 
changes is weak. On the other hand collective motions 
i.e. structural relaxations, occur through a broad distri- 
bution of very slow modes called a modes. This modes 
are very sensitive to temperature changes. In general, 
the characteristic time of the a relaxation exhibits such 
a huge increase when the temperature is lowered, that it 
seems to diverge, up to experimental evidence. For low 
enough temperatures, such that the relaxation time is 
larger than the experiment time by many orders of mag- 
nitude, the systems present some strange non-equilibrium 
features. One of the most striking among them, is the so 
called aging phenomena. It consists in a drift of the a 
relaxation time distribution, towards longer and longer 
relaxation times. This distribution keeps on evolving for 
any experimental time scale. The simplest way to char- 
acterize this drift is to perform a quench from a high 
temperature state where an equilibrium distribution can 
be achieved to a low temperature state in the glassy re- 
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gion [|[ H]. In this case, the age of the system is de- 
fined as the time elapsed since the quench, i.e. the time 
spent in the glassy phase. In theses conditions, it has 
been proven theoretically and experimentally alike that 
any measured quantities depends explicitly on the age of 
the system ||]. More precisely, the response and corre- 
lation functions do not only depend on the elapsed time 
t since the beginning of the measurement as it would be 
the case for an equilibrium system but they also depend 
on the age of the system t^^ at which the measurement 
started. Consequently, any correlation function g writes 
g{tu,,t -\- tyj). Predictions and experimental results show 
that, despite the absence of equilibrium, there is a regime 
often called asymptotic regime, where the only relevant 
time scale is the very age of the system. As a conse- 
quence response or correlation functions can be rescaled 
with t^. One finds that g{tw^t -\- 1^) — ^(^^w~f^) where 
h is a function that depends on the system [O. When 
g{tu,^t -\- tyj) = g{j-)^ the scaling is usually called full 
aging. We point out, that the existence of such a rescal- 
ing in this asymptotic regime shows that the distribu- 
tion of relaxation times drifts in a self-similar way. Such 
an asymptotic regime can only be reached if the system 
is left at a constant temperature for a sufficiently long 
time in the glassy phase. However, if some additional 
energy is transiently provided to the aging system, the 
dependance of the relaxation functions with the age of 
the system becomes extremely intricate. For example, 
temperature ramps with stops lead to the very remark- 
able " memory effect" as observed in spin glasses - see eg 
|]l3| - and recently in polymer glasses The shape 

of the relaxation times distribution is then modified in 
a non-trivial way during potentially very long transient 
regimes. The analysis of this behavior would help to un- 
derstand the precise modification of the relaxation time 
distribution by external parameters, and thus to get an 
insight of the system's internal dynamics . This is the 
scope of this paper. We have chosen to work on a col- 
loidal glass where it is easy to measure a time-dependant 
correlation function. 

For colloids, temperature is not a practical parameter. 
But since it is believed that temperature and shear may 
act similarly in theses systems we studied the in- 
fiuence of the shear upon the microscopic dynamics of a 
colloidal suspension. 

In this paper we will focus on the modification of the 
shape of the distribution of a relaxation times Pijoc). 
We will show that injecting transiently some mechani- 
cal energy into the system modifies the distribution of 
relaxation times in a dual fashion. Conversely to what 
could be intuitively expected, it can both rejuvenate or 
overage the system. In practice, we study a dense sus- 
pension of purely repulsive polystyrene beads both elec- 
trostatically and sterically stabilized. We record the fiuc- 
tuations of a laser light, multiply scattered by the sam- 
ple. We compute the two times intensity autocorrelation 
function g2{t -\-tw,tw)' The variation of g2 is directly re- 
lated to the mean square displacement of the scatterers. 



In order to "quench" the system properly, the sample 
is first presheared by a large oscillatory strain. When 
the shear is stopped we record g2{t -{-t^^t^) for evenly 
spaced t^. We show that it displays the properties of 
full aging. Then the sample is submitted to a burst os- 
cillatory strain of different amplitude and frequency. We 
monitor the change in the shape of g2 and interpret it in 
terms of changes in the distribution of relaxation times. 

III. SAMPLE PREPARATION AND 
EXPERIMENTAL SETUP 

The sample is a commercial suspension of polystyrene 
spherical beads of diameter 162 nm copolymerized with 
acrylic acid (1%) that creates a charged corona stabiliz- 
ing the microspheres. The corona prevents both from ag- 
gregation by steric and electrostatic repulsions and from 
crystallization. The suspension was carefully dialyzed to 
a polymer volume fraction of ^ = 50% ± 0.5%. The vol- 
ume fraction was determined by drying and was chosen 
to fulfill the two following criteria: 

- The sample has to be concentrated enough so that both 

a and P modes exist. The amplitude of the f3 mode 
must be sufficiently small so that the alpha mode 
can be accurately measured. The a modes must 
also display some aging behavior over all the ex- 
periment time scale. 

- The volume fraction has to be low enough so that 

the strain is homogeneous, at least on macroscopic 
length scale all through the sample. We checked 
that no effect like shear banding occurs for the con- 
centrations in use. 

Because of the divergence of the viscosity with the volume 
fraction in the vicinity of the glass transition, the suitable 
range of concentration fulfilling these two requirements 
is narrow (between 50% and 52%). 

Fig shows the experimental setup. The sample is il- 
luminated using polarized light from an argon-ion laser 
operating at a wavelength of A = 514 nm. The laser 
beam is expanded to a diameter of approximately 1 cm 
and is incident on the sample cell. Multiply scattered 
light is collected by a 50-mm Nikon camera lens and 
split into two parts. The first part is focused onto an 
iris diaphragm. The lens is set up so that a one-to-one 
image of the scattered light from the output plane of 
the sample cell is formed on the diaphragm. The scat- 
tered light emerging from the diaphragm is detected by a 
Dalsa CCD camera, model CAD1-128A, which is placed 
approximately d 15 cm behind the diaphragm. It is 
an 8-bit camera which we run at 500 frames per sec- 
ond. The images are transferred to a computer running 
at 500 MHz using a National Instruments data acqui- 
sition board, model PCI1422. The combined speeds of 
the computer and data acquisition boards are sufficiently 
high to allow data to be analyzed in real time. The ac- 
cessible times range from 10~^s to 10^ s. 
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FIG. 1: Experimental set-up. The sample is placed in a shear 
cell. The emerging light is simultaneously analyzed by the fast 
correlator and the camera. The problem of ergodicity is solved 
by performing a spatial average for the camera part, and mov- 
ing a frit glass in front of the PMT detector for the fast de- 
tection. 



The second part of the light is shone onto a moving frit. 
Its motion is controlled by a piezo actuator. The light 
emerging from the frit is detected by an optical fiber, am- 
plified by an ALV photomultiplier tube (PMT) and ana- 
lyzed by a Flex correlator. The moving frit has the same 
effect as the second cell in the two cells technique 
but we found it easier to implement. It allows an exact 
determination of the non ergodic correlation function at 
very short times (10~^5, 10~^s). A complete description 
of the technique can be found in Q . In the transmitted 
geometry the correlation functions calculated with the 
camera and with the PMT can be overlaid [|l^, |l^ . The 
dynamics of the beads can thus be probed over 12 decades 
in times. The sample is placed in a custom-made shear 
cell consisting in two parallel glass plates with a variable 
gap. For all presented experiments, the gap was set to 
1.3 mm. Oscillatory straining was realized by moving 
the bottom plate thanks to a piezoelectric device. Shear 
strain from 30% to 0.04% could be possibly applied at 
different frequencies ranging from 0.01 Hz to lOHz. 

IV. EXPERIMENTS 

A. rheological quench 

In order to check if an oscillatory strain is able to en- 
tirely rejuvenate the sample, we submitted the suspen- 
sion to a series of shear strain of different amplitudes at a 
fixed frequency of IHz for 100s. The measurement of the 
correlation function starts when the strain is stopped (see 
fig ^). The shear cessation is taken as the origin for the 
age tw of the system. For a strain amplitude above 20% 
the correlation functions become insensitive to the strain 



amplitude and to the past history. A reproducible state 
is reached. We have checked that this state does not de- 
pend on the duration of the strain application provided 
that it is superior to 40s. Typical curves are plotted on 
fig 0b for a volume fraction of 51%. 
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FIG. 2: Normalized intensity autocorrelation functions. The 
fast decrease of the function corresponds to the f3 modes and is 
measured via the correlator. The long time decay corresponds 
to the a modes -structural rearrangements. It never reaches 
any steady state on our experiment time scale. The f3 modes 
are insensitive to shear and are in a steady state. For sake of 
clarity only one curve was plotted for the short times. Inset: 
versus tnj • We find that the relaxation time varies linearly 
with the waiting time. 



The correlation function shows a characteristic two 
steps decay as earlier mentioned. A first decrease hap- 
pens around Tf^ = 10~^s. It is insensitive to shear (for 
sake of clarity only one curve has been plotted). Then 
the correlation function plateause at a value gpiat — 0.5. 
We use a regular algorithm to extract the mean square 
displacement of the beads for this plateau value. It corre- 
sponds fluctuation of position over a distance of ^ 9nm. 
Finally a second decrease takes place at a time that 
varies with the age of the system, is arbitrarily defined 
so that g{tw^tw + Tq,) = \gpiat' The inset of fig || shows 
that Toe evolves proportionally with the sample's age t^. 
All the aging part of the curve can be rescaled by t/ty^ 
as displayed on fig ||. The aging part displays thus the 
characteristic scaling of full aging. It means that Pija) 
evolves in a self-similar manner. The high shear pro- 
vided the system enough energy to entirely rejuvenate it. 
It thus corroborates the macroscopic observations previ- 
ously mentioned. It also emphasize the analogy between 
temperature and strain as far as quenches are concerned. 



B. Effect of a moderate oscillatory shear 

We now examine the influence of the strain amplitude 
upon the microscopic dynamics. All the following curves 



FIG. 3: a relaxation of the intensity autocorrelation function 
vs The rescaling of the curves is satisfactory. 



are taken in backscattering geometry for a volume frac- 
tion (j) = 50% and a measured value of 5 6Snm. As 
demonstrated in the previous section, a high shear is able 
to entirely rejuvenate our system. Hence we first submit 
the sample to an oscillatory strain of 30% for 100s in or- 
der to obtain reproducible results. Secondly, the system 
is left at rest for 10s then it is submitted to a second burst 
of oscillatory strain of different amplitude 70 and dura- 
tion du. We set the frequency to IHz. The strain history 
is displayed on fig We now take the origin of the sys- 
tem's age just after the second burst. This convention 
is purely arbitrary but it allows an easier representation 
of our results on a log scale. We record ra(7o,^w) as a 
function of for different amplitudes of 70. We take 
the curve for 70 = 0% as the reference curve. On fig 3b 
we plotted the ratio R = ra{'yo,tw)/Ta{0,tw). If < 1, 
then the effect of shear is to rejuvenate the system, that 
is to say that the sheared sample has a quicker dynamics 
than the unperturbed one. If > 1, then the inter- 
nal dynamics is slowed down by the shear application. 
We call this situation " overaging" . If a single oscillation 
is applied at IHz for Is, the sample is partially rejuve- 
nated whatever the strain amplitude may be. This is 
shown on fig ^ where all the curves lie in the region 
where R < 1. One can wonder if the partial rejuvena- 
tion corresponds to a simple backward shift of P{ra)- In 
other words, does the rejuvenation process under shear 
is the time reversal process of aging? If this would be 
the case, one could define an effective age tg// such that 

Ta(70,^w +teff) = Ta{0, tu,). 

Experimentally we could not find such an effective age 
and time translation proved inefficient to collapse all the 
^ailOitw) onto a master curve. It thus mean that P(r) 
is not modified by the strain in a self-similar manner. 
However, at long all the curves merges to R = 1, indi- 
cating that at long times the self-similarity is recovered 
as expected. 

However if 100 oscillations are applied for 100s at IHz, 
the situation is changed as shown on fig ^c. For the 
highest shear amplitude the ratio R remains constantly 
below 1 consistently with the results of the previous sec- 
tion. However for the smallest strain amplitudes, R is 
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FIG. 4: a) Strain history 7(t). We vary both the amplitude 
and duration of the second burst, b) Normalized relaxation 
time ra{jo,tu,)/ra{0,tw) vs til) after a burst of duration Is at 
IHz for different amplitudes 70 = 2.9%(«),7o = 5.9%(a),7o = 
7.9%(n),7o = 11.7%(x),7o = 14.5%(o), and complete rejuve- 
nation (■). c) Same curves for a duration of 100s. Notice that 
for the lowest shear amplitude overaging occurs. 



first inferior to one but then becomes superior as in- 
creases. The effect of shear for this amplitude is thus 
dual. Shortly after the burst the dynamics is acceler- 
ated. But after a while it becomes slower than that of 
the unperturbed case. It thus means that the transient 
shear application has modified the distribution of the re- 
laxation times not only in a non self-similar manner but 
also in a non monotonic way. 

Notice, however that at very long times all the curves 
seem to converge to R=l. It shows that regular aging 
is recovered. By considering only Tq, in the relaxation 
process we have reduced the relaxation to a single time. 
However, the shape of the correlation function is the re- 
sult of the whole distribution of relaxation times. It is 
thus interesting to compare the full correlation function 
^2(^5^^,70) to the reference one ^2(^7 ^w,0). 

Fig H shows the reference curves 5'2(^,^w,0) for tuj = 
0.1 5 (dark line) and = 60s (grey line). We em- 
phasize the fact that is now referenced from the 
cessation of the burst. The symbol curves represent 
g2{t,tw,7.9%) for = O.I5 (o), = Is (a), and 

= 60s (□). The comparison of the curves for = 0.1s 
reveals that ^2(^,^^,7.9%) starts decreasing earlier than 
92{t,tuj^0%). It is thus an indication that faster re- 
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FIG. 5: Correlation functions obtained after the second burst 
in the reference case for t^. = 10"^s (Bold line) t^. = 60s 
(Grey Line) and in the case 70 = 5.9% for tw = 10""^ s (O), 
tw = Is (A), t-w = 60s (□). Notice the crossing of the curves 
occurring around 30s. The change in the shape of the corre- 
lation function is clearly visible. 



FIG. 6: Normalized a relaxation time for 70 = 7.9% at various 
frequencies: O.lHz(o), lHz(+), 5Hz (□), lOHz (a). Notice 
that the importance of the overaging regime increases with 
the frequency. 



D. experimental conclusion 



laxation process occurs in the sample. However, at 
long times, ^2(^,^^,7.9%) crosses and then lies above 
5^2 (^, ^w, 0%). It thus mean that at long times the re- 
laxation processes in the sheared sample are slower than 
in the reference one. This is confirmed by the relative 
position of ^2(^,^^,7.9%) and ^2(^,^k;,0%) for = 60s. 
^2(^, ^w, 7.9%) decreases more slowly than 5'2(^, ^w, 0%). 
Thus, the fast relaxation times have aged and the dy- 
namics is dominated by the slower ones. The change 
in the shape of the correlation function is the very sign 
that the relaxation time distribution has not been simply 
shifted backwards in time by a constant amount, but that 
its shape has been modified with an addition of short and 
long relaxation times. 



In conclusion, we have shown that a high shear is able 
to generate in the system a distribution of relaxation time 
P{Ta) independent of the sample history and of shear 
amplitude. We then demonstrated that a moderate os- 
cillatory strain both partially rejuvenate and overage the 
system. We deduced from that point that the rejuvena- 
tion process is not a simple backwards time translation 
of P{ra). It involves a more sophisticated process during 
which the shape of the distribution is modified. However, 
independently of the imposed perturbation, P{ra) seems 
to converge to the same time-dependant distribution at 
very long times. In addition, if an oscillatory strain is 
applied similarly but at a higher frequency the overaging 
process is amplified. 



V. SGR MODEL 



C. effect of frequency 

Up to now, we have only focused on the amplitude 
and duration of the second burst. However, one can also 
expect that its frequency plays a role in the rejuvenation- 
overaging process. Indeed, one could expect that not only 
strain but also strain rate is an important factor in this 
mechanism. We set the amplitude to 7.9% and the du- 
ration of the burst to 10s. We varied the frequency from 
O.lHz to lOHz. Hence the sample was submitted to a 
number of oscillationsn ranging from 1 to 100. The result 
are displayed on fig ||. Notice that the overaging effect is 
all the more pronounced that the frequency is high. We 
believe though that the overaging behavior cannot be ob- 
served under steady strain rate as usually performed in 
rheological measurements. However, we could not check 
this hypothesis since our set up does not allow continuous 
straining. 



In a previous paper we emphasized the similar- 
ity of this results with the predictions of the simple trap 
model solved with a step in temperature. This qual- 
itative agreement reinforces the similarity of shear and 
temperature on the particles level. We also pointed out 
that a resolution of the Sherington-Kirkpatrick model for 
spin glasses in the transient regime following a temper- 
ature step gives the same results |l^, However, the 
equivalence between shear and temperature increase re- 
mains an hypothesis. We present here some results on a 
model where the macroscopic shear is coupled to the mi- 
croscopic. There are only two of such models that we are 
aware of. One is based on a mode coupling approach and 
has only be solved in the asymptotic regime |2^. The 
other one is inspired from the trap model and is called 
the Soft Glassy Rheology (SGR) model @, |, §. We 
will now solve this model for the transient strain history 
applied to our system. 

The basis of the model are the following: the system 
is described by fictive independent particles moving in a 
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fixed energy landscape. Only local minima are consid- 
ered. The particles are trapped in wells of depth E from 
which they escape in a "thermal" like fashion. The es- 
caping probability is proportional to exp[— where x 
plays the role of the thermal energy. The macroscopic 
strain is introduced as an external field that shifts the 
minimum energy levels E by an amount —^kP{t) where 
l{t) is the local accumulated strain. Thus straining helps 
to hop outside the wells. In order to calculate l{t) the 
following assumptions are made: 

i. Each time a particle escape from a trap, it falls 
in an unconstrained state of depth E with / = 0. 
The probability of falling in a trap of depth E is 
proportional to the density p{E) of trap E. 

ii. The strain rate is homogeneous all over the sample. 



iii. The two preceding hypothesis allow to define l{t) 
as the integrated local strain since the last hoping 
event U for the particle: l{t) = J^^ jdt. 

iv. The local elastic modulus k is independent of the 
trap depth. 

V. p{E) is exponential: p{E) = ^ ex.-p[—E/xg]. 

A complete description and justification of the model can 
be found in 0. 

We call P{EJ^t)^ the time-dependant probability for a 
particle to be in a well of depth E with a strain 1. The 
evolution of P{E^ l^t) in the SGR model reads as: 



^liElMl = _ p^E, I, i)e-(^-4'='^(*))/^ + rmmi) - 7^^^ (1) 
at V ^ / V ^ / al 

escaping term entering term 

external advection 



where the time unit has been set to 1, and with 

/OO nOO 
/ P{E' ,1' ,t)e-^^-"^^''"^'^^^^dE'dl' 
-OO io 

Each {E^ I) state is thus associated with a relaxation time 
te oc exp[(£^ — l/2kl'^) / x\. In absence of shear, the local 
strain remains always equal to zero, and the model re- 
duces to the Bouchaud's trap model: for x > x^, a steady 
distribution for P{E) exists. Oppositely, for x < Xg no 
steady distribution exist and P{E) reaches an asymp- 
totic regime with a full aging behavior. Subsequently Xg 
is assimilated to a glass transition temperature. If a con- 
stant shear rate is applied, the system always displays an 
equilibrium distribution of relaxation times. The model 
describes qualitativelv well many rheological features of 
soft glassy systems [0, ^ • 

We numerically solve the equation |l| by discretizing the 
equation into 100 energy levels and 100 strain level. We 
checked that the obtained results do not depend on the 
numbers of levels we use. We take k = 2^ Xg = 1 and 
X = 0.5. We remark that the values of 7 can not be 
compared to that really used in the experiments. We 
use the initial conditions of a deep quench. The sys- 
tem is supposed to have its equilibrium distribution for 
T = 00. Hence, P{E,l,0) = p{E)d{l). The external 
strain is taken as follows: 



f if t < tatt 
lit) = { fit) if tatt < t < tatt + du (2) 

[ if t > tatt + du 



f(t) is a function of t that will be specified cases by cases. 
We discuss the effect of shear by computing P(£, t) with 
e = E—^kP. Notice that the relaxation time distribution 
P(log(r£;)) (X Pis). Hence Pis) and Pirs) have the same 
physical meaning. For sake of clarity, we will use P{e) 
for which all the described effects are more visible. 



A. one step 

Firstly, we exemplify the effect of a single square pulse: 
fit) = 70. 

We first let the initial distribution P{E,l,0) = piE)S{l) 
evolves between t = and t = tatt- We chose tatt so 
that P(£^, tatt) has nearly its self-similar shape. In this 
regime, the t/t^ scaling leads to the following property 
for the distribution function: 

P{E-e,l,t) = P{EJ,t-teff) (3) 

with teff = t{l — exp[— e/x]) and e is any positive con- 
stant. For t = tatti the effect of the step is to shift the 
states (£^,0) to the states (£^,70). 

For tatt < t < tatt + du the population is now split into 
two 1 levels as the particles always jump from a (£^,70) 
state to a {E^ 0) one. One can show that, because of equa- 
tion H, the evolution of P{s, t) during the step is the same 
than that of the unperturbed case but shifted in time by 
teff = tatti^ — ^/c7q]). However, according to equa- 

tion Q the states (^,70) do not age but simply becomes 
depopulated. Only states {E^ 0) have a non zero entering 
rate. Hence there is no possible self-similar regime for 
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P(E,l,t) until P(£',7o) = 0. Hence at t = tatt + du the 
effect of shear is to shift backwards all the levels: {E^ 70) 
states are shifted to (E^O) and (E,0) states are shifted 
to — 7o) states. However, the effect of the strain shift 
differs from a simple shift in time because P{E^l^t) is 
not in its self-similar regime - despite the fact that P{€) 
has its self similar shape-. 
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FIG. 7: Distribution function P{e) when no shear is applied 
(stray line) and when a strain square pulse (70 =) is applied. 
The inset is a zoom on the region where the overaging is 
visible. 

Fig l^hows that it leads to split of the P{s) into two 
bumps. Low energies are mainly populated by the par- 
ticles that hoped during the square pulse. Conversely, 
high energies are mainly populated by particles that did 
not hop. But in addition, and as a result of the whole 
history these levels become overpopulated by jumps from 
the low energy levels. Finally, a splitting of the distribu- 
tion results from the acceleration of the kinetics of low 
energy levels and overaging originates from this splitting 
of the energy population. 

Discretizing the strain history as very short steps, one 
can solve equation ^ for any strain history. We will see 
in the following section that an oscillatory strain ampli- 
fies deeply the splitting of P{e). 

B. sinusoidal strain 

1 . effect of amplitude 

In this section we take f{t) = '^Qsin^uot) with uj = |7r. 
It corresponds to a burst of 10 cycles. Fig || shows P{e) 
0.1 tu after the burst. 

The solid line corresponds to the unperturbed case: 
70 = 0. The curves lies in increasing order of strain am- 
plitude 70 from bottom to top on the right hand-side 
of the figure. It appears that for large amplitudes -eg 
7o = 4.5- the system is accelerated. Indeed, the distribu- 
tion P{e^ 4.5) lies over the reference one for small e i.e. 
small relaxation times, and below the reference one for 
large e i.e. long relaxation times. It thus corresponds 
to a case where the system was rejuvenated. Interest- 
ingly, for moderate amplitudes. The distribution splits 




1 2 3 4 5 



FIG. 8: Calculated distribution P(£) at 5Hz for different 
strain amplitude 70 = (Bold Line), 70 = 1.7 (o), 70 = 2.5 
(a), 70 = 4.5 (□). Notice the splitting of the distribution and 
the surpopulation of both the short and large e compare to 
the reference case 70 = 0%. 



into two bumps. Both low and high e lie above the ref- 
erence curves, whereas the population of moderate e are 
depleted. In order to better understand how the split- 
ting of the correlation function influences the dynamics, 
we calculated: 

n pOO / 1 /2 

C{t^^t,t^) = J J P{E\r,t^)e^-'-'~^^~''' ^^^^E'dl' 

C(t ^ tw^tw) is the probability that a particle has not 
change its position between time and + t- It 
is actually similar to our measured correlation function 
5'2(^+^w, ^w).More precisely, ^2 is a monotonic function of 
C, going from 1 to 0, while C is going from 1 to 0. The 
detailed shape of g2{C) depends on the set-up optical 
characteristics and is of no interest in the following. Op- 
tically the situation is similar to the one of structurally 
evolving foams and is discussed in |^ . 




0.01 0.1 1 10 100 



t /tu 

FIG. 9: Same curve that figH but for the calculated C(tw^t + 
tw) for tio=0.01tu (o),6 tu (a), 50 tu (□). The amplitude of 
the strain is 70 = 2.5. 

Fig ^ displays C{t ^ twitw) for 70 = 2.5 and tw = 
0.01,6,50 tu after the burst. The reference curves for 
the unperturbed case at tuj = 0.01 -resp 50 tu- are plot- 
ted in dark -resp grey- solid lines. For 70 = 2.5, the 
correlation function C{t + ty^^ty^) for = O.Olt-u (o) 
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starts to decrease before the unperturbed curve. This is 
a consequence of the overpopulation of low energy states. 
However, the rate of decor relation is smaller because the 
depletion of intermediate energy state. Since the high en- 
ergies are overpopulated, the two curves crosses. At long 
time, the unperturbed curve lies under the perturbed 
one. For longer (□), the low energy states have aged 
and the shape of the correlation function is dominated 
by the population of the high energy. The correlation 
function for 70 = 0.06 looks 'older' than the reference 
one. Notice the excellent qualitative agreement between 
fig I and fig|. 



2. effect of frequency 

In this section we now focus on the effect of the burst 
frequency upon the dynamics. We keep the same initial 
conditions but have co varying from 27r/50 to 27r. It cor- 
responds to a number of oscillations ranging from 1 to 
50. As explained in the previous section we calculated 
C{t -\- t^^t^) for different after the burst. 




FIG. 10: Calculated To.5{tw,'jo)/ro.5{tw,0) for different fre- 
quencies of oscillation. The amplitude 70 is 10. The frequen- 
cies are O.lHz, 0.2Hz, 0.3Hz, 0.4Hz, 0.5Hz from bottom to 
top. Notice that the overaging effect is more sensitive to the 
frequency for the model than for the experiment. 



Following the analysis procedure we used in the exper- 
imental section, we define To.5(t^(;, 70) such that C(ro.5 + 
tw^tw) = 0.5. Figure^ displays To.5(t^(;, 70) as a function 
of tw ' One can see that it qualitatively predicts the same 
frequency behavior than that displayed by our sample - 
see fi^. Actually oscillatory shear strain, by mixing in 
strain and thus in £ the population, considerably accel- 
erates the dynamics of the small energy levels, while it 
does not affect the high energy levels. Thus increasing 
the frequency f{t) - for the same duration of application 
- mainly accelerates the low energy dynamics, by mak- 
ing the strain mixing more efficient. Thus the oscillatory 
strain is qualitatively similar to the square pulse strain. 
But the amplitude of the splitting in the e distribution 
and thus of the overaging considerably increases with the 
frequency. 



C. discussion 

We want first to emphasize the following point: overag- 
ing is a notion associated with a transient perturbation 
and a comparison with a reference case. It requires not 
only an analysis of the average relaxation time but of the 
global shape of P{e^t). Overaging comes from the over- 
population of the long relaxation times after the shear 
application. However, the experiments and the model 
show that the short relaxation times are simultaneously 
overpopulated. This is what we called rejuvenation. Be- 
cause of the simultaneous occurrence of theses two phe- 
nomena on different time scale, an exact analysis can only 
be performed if the complete distribution P{e) is stud- 
ied. However, the experimental time window is limited 
and the system can thus appear rejuvenated or over aged 
depending on the relaxation times that we can probe. 
We now wonder if a regime where rejuvenation alone can 
be achieved for any time scale after a perturbation 70 
and a duration du. The experiment duration are too 
limited to give a satisfactory answer. However we can 
make the following remarks about the SGR model. In 
this model, aging occurs because for x < Xg the escap- 
ing rate Tout{E) is inferior to the entering rate Tin{E) 
for high enough energy wells. ToutiE) is proportional to 
ex.-p[—E/x] and thus depends only on the on the well's 
depth. Oppositely, Tin{E) is proportional to T{t) and 
thus depends on the whole distribution P{E). Let us 
first examine the population of the level (E,0) before the 
shear such that 

E:^ Ec = ^kjQ + X ln{du) 

. When a shear is applied, ToutiE) is simply multiplied 
by a factor ex.-p[^kP /x] during a time du. Hence es- 
caping probability of this population will remain nearly 
zero despite the shear application. Theses particles will 
remain in the states {E^ 0) after the shear has been ap- 
plied. However, for these energies, the entering rate is in- 
creased because of the overpopulation of small energies. 
Overaging comes thus from the fact that the escaping 
rate remains unchanged whereas the entering rate is in- 
creased for large energies. Hence overaging in the SGR 
model occurs for any finite 70 and any finite duration. In 
this case a complete rejuvenation is impossible. However, 
the assumption that the trap stiffness is independent of 
the well depth is only a first step assumption. One could 
imagine that the deeper the well the stiffer its elastic 
modulus. One could for instance assume that: 

k = ko kE 

where z^: is a constant of proportionality. Making this as- 
sumption first changes the instantaneous elastic modulus 
< k > of the model. It is no longer ko but now 

< k >= ko ^ < E > 

For X > Xg Si stationary distribution of P(E) exists and 
we find < k >= kn -\- k . In the limit of infinite 
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temperature < k > tends towards ko + K,Xg. However, 
when X < Xg^ the average energy < £^ > is proportional 
to log{tuj) in the asymptotic regime. The instantaneous 
elastic modulus has thus a logarithmic dependance with 
the ag e of the system as observed in many systems ||, 
|2^. An other change affects the escaping rate. Tout 

is now proportional to exp[— [ '^"^ E — ^^o^^]]- The 
strain now plays a role similar to a temperature change. 
In this case, the value of Ec now reads: 

|/c7o + x\n{du) 

1 2 

One can thus define a critical strain 70c = f for which Ec 
is infinite. Hence a complete rejuvenation ca be achieved 
in a finite amount of time. From an experimental point of 
view, the distinction between the two approximations by 
a rejuvenation experiment are not yet concluding. Some 
other experiments are being performed and the conse- 
quences of this modification for the stiffness will be care- 
fully discussed elsewhere. 



VI. CONCLUSION 

In conclusion, we have shown that an oscillatory strain 
can act on a glassy colloidal suspension in a dual fash- 
ion. It can rejuvenate it by erasing all its past memory. 
It can also overage it by accelerating the aging process. 
Both processes comes from an acceleration of the rear- 
rangement rates during a transient time. The dominat- 
ing effect depends on the amplitude and duration of the 
shear. We showed that the effect is qualitatively well ex- 
plained by the SGR model. This model was solved for a 
realistic strain history. We showed that the distribution 
of relaxation time was separated into two bumps, one 
corresponding to rejuvenation and the other one to over- 
aging. The influence of this two bumps on the position 
autocorrelation function was studied. We showed that 
the model and the experiments are highly comparable. 
We believe that this phenomena is not a special feature 
of colloidal systems but can be observed in polymer - as 
indicated by preliminary results performed in our lab - 
or in spin glasses. 
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